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^ . Abstract 

OO ■ We analyse the structure of one-loop self-energy graphs for the p meson in real time formulation of finite 
CN . temperature field theory. We find the discontinuities of these graphs across the unitary and the Landau cuts. 

| ! These contributions are identified with different sources of medium modification discussed in the literature. 

' We also calculate numerically the imaginary and the real parts of the self-energies and construct the spectral 

^p" 1 . function of the p meson, which are compared with an earlier determination. A significant contribution arises 

Qh' from the unitary cut of the ttco loop, that was ignored so far in the literature. 

1 Introduction 

t^j- ' The in-medium propagation of vector mesons, particularly the p, has been extensively studied in the literature, 
as reviewed in Refs.[TJ[2]. The reason is, of course, that it controls the rates of dileptons and photons emitted 
from the hot and dense matter, created during heavy ion collisions. The recent precision measurement [3] of 
the in-medium p spectral function encourages further theoretical investigation. 

■ The effect of medium on the vacuum propagation of the vector meson is generally believed to arise from 
s * two sources pQ. One is the change in its pion cloud, given essentially by the tttt self-energy loop [I]. The other 
^p. \ is the collisions suffered by the vector meson with particles in the medium. The latter effect can be calculated, 

■ broadly speaking, in two different ways. Kinetic theory expresses this collision rate in terms of the spin average 
^ . of the squared scattering amplitudes [5]. This rate is simply related to the imaginary part of the self-energy 

^ | of the meson [6]. The effect of collisions can also be obtained from the self-energy tensor given by the virial 
formula, which relates it directly to the scattering amplitude itself [H [91 [101 [11] . In the absence of scattering 
data, these amplitudes are generally calculated from the relevant Feynman graphs. 

As shown by Weldon and others [61 I12j . the two sources modifying the free propagator find a unified 
description in terms of contributions from branch cuts of the self-energy loop. In addition to the unitary cut, 
present already in the vacuum amplitude, the thermal amplitude generates a new, so-called Landau cut. While 
the in-medium modification by the cloud of virtual particles (mostly pions) is obtained from the unitary cut, 
the effect of collisions with surrounding particles is given by the Landau cut. Thus the two sources of medium 
modification are automatically included in the calculation, if we retain the contribution of both the cuts. The 
relative importance of these cuts from different graphs depend on their thresholds, besides the couplings at the 
vertices of the graphs. 

In this work we take the one-loop self-energy graphs for the p meson and find all the discontinuities associated 
with the branch cuts. The loops are formed with one internal pion line and another which may be the pion 
itself or any of the heavy particles, namely u, a\ and hi, up to a mass of about 1.25 GeV. The resonances are 
treated in the narrow width approximation. The vertices of the graphs are obtained from chiral perturbation 
theory. The calculations are carried out in the real time version of thermal field theory. 
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In Sect. 2 we start with the two-point function at finite temperature of the vector current having the 
quantum numbers of p. Here we also review briefly the methods of chiral perturbation theory to obtain the 
different interaction Lagrangians needed to evaluate the self-energy loops. In Sect. 3 we describe the kinematic 
decomposition of the tensor amplitudes. In Sect. 4 we write explicitly the self-energies from different loops and 
separate these analytically into their real and imaginary parts. The cut structure of the self-energy function and 
the discontinuities across the cuts are obtained in Sect. 5. In Sect. 6 we evaluate numerically the imaginary and 
the real parts of the self-energy and construct the spectral function of the p meson. Finally Sect. 7 discusses 
the assumptions and the approximations entering these evaluations. Here we also compare our results with an 
earlier determination. The Appendix gives a summary of the real time theory, needed in the present work. 

2 Preliminaries 

To study the p meson propagator, we do not start directly with the two-point function of the p meson field, but 
consider instead that of the vector current, having the quantum numbers of the p meson. In the two-flavour 
QCD theory, this current is given by 

V;{x) = q{x) 1 J-q{x), ?=(d)> i = l,2,3 (2.1) 

Conceptually we then keep contact with the fundamental theory and deal with a conserved current in the limit 
of SU (2) flavour symmetry. At the same time we can address directly the physical processes, such as dilepton 
production in heavy ion collisions. 

Here we work in the real time formulation of the thermal field theory, where a two-point function assumes 
the form of a 2 x 2 matrix. Accordingly we have the matrix two-point function 

T l ^ ab (E,q)=i J d 3 xdTe iET - i ^{T c V*(x)Vj(0)) ab (2.2) 

where (O) denotes the ensemble average of an operator O, 

(O) = Tr(e- pH 0)/Tre- pH (2.3) 

and Tr denotes trace over a complete set of states. The superscripts a,b (a, b = 1, 2) are thermal indices and 
T c denotes time ordering with respect to a contour in the plane of the complex time variable r, to be specified 
in the Appendix. However, as we demonstrate there, the basic quantity is again the vacuum two-point function 

1%(q) = i f d 4 xe^(0|Ty;(x)^(0)|0) , (2.4) 

from which one may easily construct the thermal components. 

In the region of low q 2 , the p meson pole and different low mass branch points contribute to the two-point 
function, which may be calculated in terms of a few hadronic states. As q 2 rises, too many hadronic states 
start contributing and a continuum sets in. Here the QCD coupling parameter becomes small, allowing a 
quark-gluon perturbative calculation of the two-point function. 

In this work we are interested in finding the leading temperature effect modifying the free propagation of 
the virtual p meson at q 2 below the continuum. At low temperature, the medium is populated mostly by pions. 
Thus we consider one- loop self-energy graphs of Fig. 1(a), consisting of a pion and another low mass hadron h, 
along with the pion itself. There is a series of such hadronic states with increasing masses. Of course, because 
of the presence of thermal distribution function, their contributions fall off exponentially with their masses at 
sufficiently high q 2 . But it is not a priori clear, up to what mass the hadrons should be retained in our low 
q 2 calculation. Here we take the particles h to consist of ui,hi and cti, postponing discussion on this point 
to the last Section. Along with these graphs, we also include the seagull graph of Fig. 1(b), where the loop 
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Figure 1: One- loop graphs for the two-point function contributing to self-energy of p meson. 



represents a single pion propagator, those of heavy particles contributing insignificantly at low temperatures 
that we consider here. 

The two-point function in the hadronic phase can be evaluated in chiral perturbation theory by using the 
method of external fields [13) . Here one introduces (classical) external vector field v^x) coupled to the vector 
current V*(x) (and an axial- vector field a* (re) coupled to the axial- vector current A 1 (x)). The resulting terms 
perturb the original QCD Lagrangian as 

Cqcd -> Cqcd + viYt + «X ( 2 -5) 

The global chiral symmetry group SU(2)r x SU(2)l of massless QCD is thus raised to the corresponding local 
gauge invariance. To arrive at the resulting effective field theory, we have to define the various fields and their 
covariant derivatives transforming properly under the symmetry group [14} \TE[ [TB] [T7] . The Goldstone fields 
representing the pions may be described by the unitary matrix, 

3 

U(x) = e**W F * , tt(x) = Y J ^ l {x)T i , (2.6) 

i=i 

where the constant F n can be identified with the pion decay constant, F n = 93 MeV and r a are the Pauli 
matrices. One also defines u(x), such that u 2 = U, because of its transformation rule linking SU{2)^l to the 
unbroken group SU(2)y- The non-Goldstone iso-triplet field R fl (x) representing p and a± mesons, 

1 3 

R ll {x) = ^=Y. R » T ^ ( 2 ' 7 ) 
V2 i=i 

transforms as 

R^^hR^ h\ /i(vr) G SU(2) V , (2.8) 

while the singlet field S^x) representing u and h± mesons, of course, remains unchanged. The covariant 
derivatives of U and R^ are defined by 



and 

with the connection 



DpU = d^U - i( Vfl + a^U + iU(v„ - afl ), (2.9) 

V lt R u = d tL R l/ + \r„R u ] (2.10) 

t p = \ (yidv - K v n + a p)\ u + u \ d P ~ K v » ~ a p)W ) (2- 11 ) 

Finally we have the usual definition of field strengths F^ L associated with the non-abelian external fields 

(v/j, ± ci/j), where = ^ u u r V2 an d o-fi = ]C a /i r V2- We can now write the leading term in the effective 
lagrangian for the Goldstone bosons as 

C=?f[(D l tUD»U)+ml(U + U j <)} (2.12) 
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It is possible to define new variables related to D^U and F^ L , namely [IHl [T7] . 

u p = iu^DpUv) = ul (2.13) 

and 

/£" = ±v*Fg'u + uF£ v v) (2.14) 

which transform exactly as Rp(x) in (2.8). We now have all the variables relevant to non-Goldstone bosons, 
that transform under SU(2)y only. It is now simple to write the different pieces of the Lagrangian density, that 
are invariant under SU(2)y and hence also under SU{2)r x SU(2)l- Accordingly, we have the kinetic part as 

Ckin = £ (-{(R^n + \ml(R»Rn) + £ (- \ S,uS^ + \mlSpS^ , (2.15) 

where 

Rnv = V M i2„ — VuRfj,, Sp V = d u S v — d v S^ L . (2-16) 

and the two sums run over the iso-triplets and the iso-singlets. Here and below the symbol (A) stands for the 
trace of the 2x2 matrix A. The interaction vertices to leading order in powers of derivatives (momenta) and 
external fields are given by [HU dTJ [18] 

(2-17) 

Each of these interaction terms is typically a series in powers of the derivative of the pion field. Retaining terms 
to lowest order, we get the interaction vertices as, 



m p 



1C 

lk d »P» ■ d ^ x 

m p F4 

-^-h^(dpp u - d v pp) ■ d u TT 

+§-(d„p v - d u p u ) ■ x d v v (2.18) 

The magnitude of the coupling constants may be determined from the observed decay rates of the particles 
PI]. Thus the decay rate T(p° -> e+e") = 6.9 keV gives F p = 154 MeV. The decay rate T(p 2vr) = 153 
MeV gives G p = 69 MeV. Similarly the decay rates Y(oj 3tt) = 7.6 MeV, T(hi pir) ~ 360 MeV and 
T(ai — > pir) ~ 400 MeV give respectively g\ = .87, #2 = 1-0 and 53 = 1.1. 

The leading term for the four-point vertex needed in the seagull graph (Fig. 1(b)) arises from the kinetic 
term for the p meson appearing in (2.15), 

A ( nl = -^<[/A^]k,d^]> (2.19) 

We note that this term in the chiral Lagrangian reproduces the famous current algebra result for the leading 
threshold behaviour of pion scattering off a heavy target, which is p in the present case |20| . 
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3 Kinematics 



The isospin structure of T p b ' 1 ^ is clearly given by 5 l i , which we omit from now on. We need to analyse its 
thermal and Lorentz index structures to extract the relevant scalar functions. 

We introduce the (four-dimensionally) transverse p meson propagator G°^ v by removing from T p b the factor 
K = (F p q 2 /m p ) 2 , representing the current-p coupling at both ends of graphs of Fig. 1, 

2#(g) = KGfM (3-1) 

The free elements are 

T$«\ q ) = KG$ ab (q), GjW«*( ff ) = {-g, v + ^) D°»(q), (3.2) 

where D ab (q) is the free thermal propagator of a scalar particle of mass m p . Comparing with Eq.(A.15), we 

find that G^} ah (q) is indeed the (four-dimensionally) transverse part of the free vector propagator D p b u (q). 
Now G satisfies the familiar Dyson equation, 

GfM = G^ ab (q) ~ Gf x ac (q)H^{q)G db M (3.3) 

where IT^, (g) represents the p meson self energy. As shown in the Appendix, we can diagonalise the thermal 
matrices to get rid of the thermal indices, getting the Dyson equation in terms of the diagonal elements (denoted 
by a bar) as 

GW?) = ^2(?) - G^(q)U Xa ( q )G av ( q ), G^(q) = (-g^, + ^) g2 _ ~ 2 + ie ( 3 " 4 ) 

Here IL fll ,(q) is related to the 11-component by Eq. (A. 18) of the Appendix. 

Being free from thermal indices, the Dyson equation may be solved in a way similar to that in vacuum. The 
vector current conservation leads to 

q fi G^(q) = q v G^(q) = 

At finite temperature there are two independent second rank tensors satisfying this condition. To construct 
them we introduce the four-velocity u p of the heat bath (not to be confused with the field u p defined in the 
previous section), even though we shall do explicit calculations in its rest frame [u 11 = (1,0)]. We thus have 
one more scalar variable, u • q in addition to q 2 . Introducing also the transverse variable u p = u p — u ■ qq^/q 2 , 
we can build two independent tensors 

Pfiv = -9^u + -^2 =2 V"<" Qp-y = -^p~u p u u , q = (u-q) - q , (3.5) 

satisfying the projection properties 

PP = -P, Q-Q = -q 2 Q, P-Q = (3.6) 

While both P and Q are four-dimensionally transverse, P is also 3-dimensionally transverse. In the literature 
one generally finds the factor q 2 instead of (q 2 ) 2 in the definition of Q pu . However, at finite temperature 
dynamical singularity can appear at q 2 = 0. The additional factor of q 2 keeps the kinematic covariant regular 
at that point. 

We now decompose G pv as 

G pv = P pv Gt + Q^iuGi (3-7) 

One can further show that the self energy tensor Ii pv is also transverse |21| . leading to an identical decompo- 
sition, 

IV = P pi M t + (3.8) 
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Inserting these decompositions in the Dyson equation (13.4p we get the solution, 

G t (q) = ^ Gi(q) = \- ~ l (3.9) 

q z - m z - IL t {q) <T q z - m z - q z Ili{q) 

The scalar self-energies may be obtained from the tensor one by taking trace over its indices and contracting 
them with the velocity four vector, 



n t = -i(nJ + ^n 00 ), n, = ln 00 , n Q0 = ^n M , (3.10) 

2 F q z q z 



Finally we note a kinematic relation between the transverse and the longitudinal components of the prop- 
agator at q = 0. As q — > 0, the kinematic structures Py and Qij depend on how the limit is reached. This 
unphysical dependence is eliminated, if we take 

G t (q ,q = 0)=q%G l (q ,q = 0) (3.11) 

Clearly a similar relation must also hold between H"^, which is already implied by Eqs. (3.9). 



4 Evaluation of loop integrals 



One-loop graphs were evaluated in the imaginary time formulation by Weldon [6] and subsequently in the real 
time version by Kobes and Semenoff |12j . which we follow to evaluate the graphs at hand. We begin with the 
general form of the self-energy loop in vacuum with internal lines for pion and a hadron h, 



/d k 
——N^(q, k)A n (k)A h (q - k) , (4.1) 



where A n ^{q) are vacuum propagators of scalar particles of mass m n and rrih respectively. Tensor structures 
associated with the two vertices and the vector propagator are included in Nau- Introducing the following three 
gauge invariant (transverse) tensors, 

A a p(q) = -g a p + q a qp/q 2 , 

B a p(q,k) = q 2 k a k/3 - q ■ k{q a kp + k a qp) + (q ■ k) 2 g al3 , 

C a p{q,k) = q 4 k a kp-q 2 (q-k)(q a kp + k a q/ 3 ) + (q-k) 2 q a qp. (4.2) 
we get the tensor N pu for the different graphs of Fig. 1(a) as 

Nfr\q,k) = -4 (jL V {B, v + ftf V) 
N^H^k) = -(f)\B, u --^ r C, u ) 

\ r K/ " l h 1 

Nfr){q,k) = -2 Q^V (B^ - -A-C^) (4.3) 



»i 
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At finite temperature the 11-component of the corresponding self-energy matrix is given by 



d 4 k 

-^N^{q,k)A w (k)A h (q - k) 



d 4 k 

NUq, k){A h {q - k)n(uj)5(k 2 - ml) + A n (k)n(uj')5({ q - kf - m 2 h )} 



(2vr)3 ■ 
d A k 



j^N^q, k)n(uj)n(uj')6(k 2 - ml)S((q - k) 2 - m 2 ) , (4.4) 
on inserting the expression for the thermal propagators D^ h {q) as given by fA8l) . Here n is the distribution 



function at temperature T and uj and uj' are the energies of the pion and the hadron, uj = ym 2 + k 2 , uj' 



m h + (q ~ k) 2 . Here the first term refers to vacuum. The second and the third terms are medium dependent, 
containing the distribution function linearly and quadratically respectively. The third term is purely imaginary. 

It is simple to evaluate the k$ integral of each of the three pieces in (4.4) and read off the imaginary part. 
Then we get ImTI(g) on using (A. 18), associated with the time-ordered product. The imaginary part associated 
with the retarded propagator is also of interest, particularly in writing the dispersion relation for H(q). The 
two imaginary parts are the same, except for a factor of e(qo). Keeping the retarded propagator in mind, we 
write (5.1) below for lmll(q). For the real part of IL{q), the medium dependent contribution is obtained from 
the second term in (4.4), 

R tt (\-[ d3k / n(u) ( N^jkp = u) N^(k = -co) 
HelVW-y (27r)3 | 2u [ {qo _ Uj) 2_ ujl 2 + {qo + uj)2 _ UJ ,2 

n{(J) ( NjjJJcq = q - uj') N^jkp = gp + ui') > , 
+ 2uj> V (qo - u') 2 - co 2 (?o + ^ 2 -'" 2 n ' " }i 



where the integrals are principal valued. 

The seagull graph (Fig. 1(b)) involves a contraction of two pion fields with a single derivative at the same 
space-time point that appear in the interaction Lagrangian (2.18). As a result, the thermal components of the 
resulting self-energy tensor has the form 

/d A k 
— q .kD ab (k) = 0, (4.6) 

if we recall that D ab is a function of k 2 and k^. Thus in chiral perturbation theory, the contribution of each 
graph of Fig. 1(a) is transverse by itself, without requiring the seagull graph of Fig. 1(b) to do so, the latter 
being actually zero. 



5 Analytic structure 

In writing the expression for the imaginary part of the self-energy, we convert the tensors N^ v and Ii^ u into 
scalars by combining their components as in (3.10). Also, throughout this Section, we omit for brevity, the 
subscripts t and I of the scalars, as the equations will be valid for both the transverse and the longitudinal 
components. We thus get 

/d 3 k 
(2^W X 

[N(ko = o;){ (1 + n{ui) + n{ui ))5(qo — u: — uj') — (n(uj) — n(uj'))5(qo — uj + uj )} 

+N(k = -uj){(n(uj) - n(uj'))5(q + u - J) - (1 + n{uj) + n(uj'))5(q + u + uj')}} (5.1) 
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Figure 2: Branch cuts of self-energy function in qo plane for fixed q given by nh loop. The quantities (71,2,3 denote 
the end points of cuts discussed in the text : q\ = \J (rrih + m w ) 2 + |ef| 2 , q 2 = \/ (nih — m^) 2 + \q\ 2 and q% = \q\. For 
the 7T7T loop, q2 collapses to q%, giving only one variety of Landau cut. 

We shall label the four terms above by the indices 1, 2, 3 and 4 in the sequence they appear in (5.1). The 
functions Ntj for different graphs may be obtained from (4.3). Here it suffices to note that these scalar functions 
depend on scalars built out of the four- vectors, q^, and u^, in such a way that they remain invariant under 
the simultaneous change of sign of qo and k,Q. We thus get the symmetry relations 

Imn 3 (-g) = -ImTT 2 (g) 

Imn 4 (-g) = -ImlTi(g). (5.2) 

The distribution functions present in different terms of (5.1) may be understood in terms of decay and 
recombination (inverse decay) probabilities for processes at the vertices of the loop graphs [6]. Denote n(oj) 
and n(uj'), for short, by n and n' . Then writing, for example, 1 + n + n' = (1 + n){\ + n') — nn' , the first 
term becomes proportional to the probability for the decay p — )■ tt + h with the expected statistical weight 
(1 + n)(l + n') for the stimulated emission minus the probability for the inverse decay tt + h — > p with weight 
nn' for absorption. Similarly, writing — (n — n') = n'{\ + n) — n(l + n'), the second term gives the probability 
for p + h — > tt with weight n'(l + n) minus that for ir — > p + h with weight n(l + n'). 

The regions, in which the four terms of (5.1) are non- vanishing, give rise to cuts in the self-energy function 
(Fig. 2). These regions are controlled by the respective 5- functions [22J . Thus, the first and the fourth terms 
are non- vanishing for q 2 > (nih+m^) 2 , giving the unitary cut, while the second and the third are non- vanishing 
for q 2 < (m/j — m^) 2 , giving the so-called Landau cut. The origin of these cuts is clear from the discussion in 
the last paragraph. The unitary cut arises from the states, which can communicate with the p. These states 
are, of course, the same as in vacuum, but, as we see above, the probabilities of their occurrence in medium 
are modified by the distribution functions. Since it is the distribution function for pions and not for heavy 
mesons that dominates the medium dependent probabilities, the contribution of the unitary cut is referred 
in the literature as due to modification by the pion cloud in medium pp. On the other hand, the Landau 
cut appears only in medium and arises from scattering of p with particles present there. We note that this 
contribution appears as the first term in the virial expansion of the self-energy function [TJ [10] [TT] . 

We now come back to evaluate (5.1) explicitly, giving discontinuities of the self-energy function across 
the cuts in the qo plane for fixed \q\. We consider a ixh loop, the tttt loop being a special case. Writing 
d?k = 2ix\Jtjj 2 — m 2 uduj sin 9d6, where 6 is the angle between q and k, we can readily integrate over cos 6 using 
the 5-functions. But we have to take into account the physical requirement, | cos#| < 1, which, as we shall see 
presently, reduces the a priori range (m^ to 00) of integration over lo. Below we inspect each of the terms of 
(5.1) and write down explicitly only those pieces that contribute for positive go- 
Consider the first two terms of (5.1), for which we have (go — ^) 2 = w /2 5 giving 

2|9lV w ~ m i 
Then the inequality | cos#| < 1 becomes 

q 2 {uj -uj + ){oj -oj^) < (5.4) 
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where co± are the roots of the quadratic equation for u, 



u, ± = ^{q ±\q]e(R 2 )v}, v{q 2 ) = ^l-^f. (5.5) 

For the first term in (5.1), for which q 2 > (m^ + m^) 2 , as already stated, we have R 2 > and v < 1, so 
that both uj + and cj_ have the same sign as that of go- Then this term is non-zero only for positive go with the 
integration variable u restricted to u;_ < oj < lu+. Changing wtox, uj = (R 2 /2q 2 )(qo + \q\x), we get 

Imrii(go,g) = -— — o / dxN(x){l+n(uj)+n(q -uj)} , qo > V (m h + m n ) 2 + \q\ 2 . (5.6) 
3i7rg z J_ v 

For the second term in (5.1), we split the region q 2 < [m^ — m^) 2 into two segments, namely, < q 2 < 
(m/j — m^) 2 and q 2 < with positive and negative values of g 2 , so that the inequality (5.4) may be applied 
immediately. Proceeding as before, we find that the first segment leads to a cut in the negative go region 
(— y 7 (rrih — m^) 2 + |gj 2 < go < — |g|) with u restricted to w_ < w < u + , which we do not write explicitly. The 
second segment gives a cut in a region over both positive and negative values of go with oj > cj_, getting 

R 2 f°° 

Imn 2 (g ,g) = 32 ^ 2 J dxN(x){n(uj) - n{uj - g )} , - |g| < g < |g| . (5.7) 

The last two terms of (5.1) may be analysed in the same way. Here the expressions for cos and u± remains 
the same as before, except for reversal of sign of go- The position of the cuts in the go plane are thus obtained 
from the previous ones by reflecting them at the origin. Thus the third term gives rise to two pieces of imaginary 
part. With uj = (R 2 /2q 2 )(-q + \q\x), we get 

R 2 dxN(x){n(uj) - n(q + uj)}, |g| < g < ( m h ~ m n ) 2 + \q\ 2 

I m U 3 (q ,q) = —^{ J -v fO0 (5.8) 



32vrg2 j _ (°° dxN{x){n{Zj ) _ n(g +5)}, -\q\ < g < \q\ 



The fourth term contributes entirely to negative values of go- 

The functions N t j(x) for different graphs are given by (4.3) with the trace and the 00-component of the 
three tensors calculated from (4.2), 

A» = -3, A)0 = l \ (5.9) 

BH = mW + ^, B 00 = -^(l-x 2 ) (5.10) 

CZ = q\mW-^), C 00 = ^x 2 (5.11) 



6 Numerical evaluation 

We begin with the results of numerical evaluation of the different graphs of Fig. 1(a) for the self-energy of the 
p. As usual, we retain the vacuum contribution in the imaginary parts only, assuming the real (divergent) parts 
to renormalise the p meson mass. We calculate the self-energies as a function of ^/q 2 = M at fixed values of 
the three-momentum q and temperature T. It thus suffices to calculate the self-energies in the time-like region, 
for positive values of go starting from go = |g|. Then the first part of the Landau cut (0 < go < |g|) cannot 
appear in our calculation of the imaginary parts. 

The 7T7T loop is distinguished by a large imaginary part of the self-energy, its vacuum part giving T p = 

Imll /m p = 153 MeV at M = m p . Clearly it is only the unitary cut in the time- like region that gives the 
imaginary part. The results for this loop are shown in Fig. 3. 
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Figure 3: The imaginary and the real parts of self-energy from 7rvT loop in upper and lower panel respectively. The 
longitudinal and transverse components are shown separately. 
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Figure 4: The imaginary and the real parts of self-energy from different nh loops in the upper and lower panel 
respectively. The quantities are averaged over polarisation. 
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In showing the results for other loops, we average their imaginary and real parts over the transverse and 
longitudinal components, 

U( q ) = ^( q 2 U l (q) + 2U t (q)) (6.1) 



They are shown in Fig. 4. Here it is only the second part of the Landau cut (\q\ < qo < \J (rrih — m^) 2 + \q\ 2 ), 
which contributes to the imaginary part. The only exception is the irui loop, where the unitary cut (qo > 
\J (rrih + m n ) 2 + \q) 2 ) also contributes, its threshold for other loops appearing outside the range of M plotted 
here. The ttuj loop dominates up to about M ~ 500 MeV, beyond which the nai loop takes over. The rising 
trend of the imaginary part at the upper end is due to the contribution of the unitary cut. While the imaginary 
parts add up, there is appreciable cancellation among the real parts of different loops. Fig. 5 shows separately 
the behaviour of the transverse and the longitudinal components of both the real and the imaginary parts of 
self-energy, summed over ttK loops. Though they differ at low M, they tend to converge at higher M. 

Finally we come to the transverse and the longitudinal components of the spectral function, defined by 

tm (M 2 -m2-(l,g2 )SRe n tj/ )2 + {(l, g 2 )Elm n M }2 

where the summation extends over all the loops. Again we take the average over the two components, 

ImG(q) = i(g 2 ImG z + 2ImG t ) (6.3) 

It is drawn in Fig. 6 for three different temperatures at q = 300 MeV, showing a reduction in the height of 
the peak with the rise of temperature. The pole position, defined by the zero of the real part of of the inverse 
propagator, shifts from M = 775 MeV in vacuum to M = 763 MeV at T = 150 MeV. Fig. 7 shows this average 
spectral function for different values of three-momentum \q\ at T = 150 MeV. It appears to depend little on 
the magnitude of q. 



7 Discussion 

We present a detailed analysis of the singularities of a one-loop, meson self-energy graph at finite temperature. 
It is carried out in the real time formulation of field theory in medium. The branch cuts are obtained in the 
(complex) energy plane at fixed three-momentum. The discontinuities across the cuts are obtained explicitly. 
The unitary and the Landau cuts arise from what are known in the literature respectively as modifications by 
pion cloud and in-medium scattering, both effects depending largely on the pion distribution function. The cut 
structure unifies the discussion of the two parts of the self-energy function. 

We also evaluate numerically the self-energy and the resulting spectral function of p, which may be compared 
with an earlier determination by Rapp and Gale [9]. Before we do so, however, it would be useful to outline 
the differences in the two treatments of the problem. Of the series of resonances with increasing masses, which 
can contribute in the loop calculation, we keep only cj, hi and a±, while they include three more. Here we note 
that although we focus on the p meson pole with self-energy corrections, the physical quantity is the complete 
two-point function, consisting of this pole and the continuum. As we take fewer resonances to calculate the 
self-energy, we expect the corresponding continuum to begin at a lower value of q 2 than theirs. The dilepton 
spectra emitted in the heavy ion collisions depend on the complete two-point function. It would thus be of 
much interest to work out this emission spectrum with these different evaluations of the two-point function and 
compare with experiment. 

They [9] also include vertex form factors and finite widths of resonances in their evaluation. The form 
factors go to decrease the contributions from higher momenta. So do the distribution functions for calculations 
in medium. Thus while the absence of form factors may lead to a significant difference for results in vacuum, it 
is not so for their medium dependent parts. The finite widths referred to are those of resonances inside a loop 
integral, where an integration over phase space is involved. In such cases, the results with finite and narrow 
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widths generally differ by no more than, say 10%. The case of A(1237) is an example [7j. Also approximating 
the vector and the axial-vector spectral functions themselves by sharp peaks for p and a\ in the well-known 
Weinberg sum rules |231 124] . one commits an error of only a few percent [25J. 

Finally there is some difference at the level of basic formulae. Their expression for the real part of the 
self-energy due to singularities in the qo < region differs from ours (4.5)), even though the difference is 
numerically small. Also the structure of some of the interaction Lagrangians are different, but again it would 
not affect the results significantly in the neighbourhood of M = m„. 

Keeping all these differences in mind, we may expect the results of our calculation to agree with that of ref. 
[9] to within, say 20%. Indeed, the imaginary parts agree well within this expectation in a wide region in M 
around the p meson mass. At higher values of M, however, our calculation shows a substantial increase due to 
the inclusion of the unitary cut of the itoj loop. 

To conclude, we derive the discontinuities across all the cuts of a self-energy loop of a vector meson. In 
the literature [9] the self-energy is constructed from the unitary cut of the tttt loop and the Landau cut of 
different irh loops. By contrast, we find that the unitary cut of at least some of the irh loops (here ttuj) may 
also contribute significantly. 
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Appendix : Real time field theory at finite temperature 

In this Appendix we review briefly the real time formulation of thermal field theory. We construct scalar and 
vector propagators and diagonalise them, showing that such matrices are actually given by a single analytic 
function, coinciding essentially with the corresponding result from the imaginary time formulation. But the real 
time procedure requires neither the frequency sum to evaluate the loop integrals, nor any analytic continuation 
to physical energies. 

Free scalar propagator 

First consider the (free) propagator for the scalar field (j)(x), 

D(t, t'; x-x') = i{T c 0(T, x)4>{t\ x')) . (A.l) 

where x^ = (t,x). Here r and r' are any two points on a contour in the plane of the complex time variable. 
The contour begins at — T say on the real axis and ends at — T — i/3, nowhere moving upwards and remaining 
entirely within the analyticity domain — j3 < Im(r — r') < [26]. Apart from these restrictions we keep the 
form of the contour arbitrary at this stage. The delta- and theta- functions on the contour will be denoted by 
the subscript c. 

The thermal propagator satisfies the same differential equation as the one in vacuum, 

(□ c + m 2 )D(r, r'; x-x') = 5 c (t - t')5 3 {x - x'), (A.2) 
but obeys a different, thermal boundary condition. If we write 

D(t, t'; x-x) = 9 c (t - t')D + (t, t'; x - Sf ) + 6» c (r' - t)L>~(t, r'\ x-x'), (A.3) 

and use translation in r and cyclicity of the trace in either of D^, we get the so-called KMS boundary condition 

as, 

D-(t,t';x-x') = D + {T-ip,T';x-x') (A.4) 
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To obtain the solution it is convenient to take the spatial Fourier transform, 

d 3 k 



D{t, r';x-f) 



(A.5) 



The r-dependence is then given by 

-72 



(J~2+^) D(r,T';k) = 5 c (T-r') 



k 2 + m? 



(A.6) 



It is now an elementary exercise in Green's function to solve Eq.(A.6) with the boundary condition (A. 4), 



(A.7) 



1 



-, getting the particle number distribution function in the propagator. 



where n(ui) = 

Of the variety of possible contours in the complex time plane [27], two are specially interesting, namely the 
closed one [28] and the symmetrical one [29J. We now choose the latter contour (Fig. 8) with T — > oo, when 
D(t — t' ,k) reduces effectively to four components, which may be assembled in the form of a 2 x 2 matrix. The 
Fourier transform D ab (k) of these components can be defined with respect to real time to give [29 1 130 1 f3T | I32j. 



ii 



■(D 



22\* 



A(k) + 2Trin5(k 2 - m 2 



D 

D 12 = D 21 = 2myjn{\ + n)5{k 2 - m 2 ) 



where A (A;) is the Feynman propagator in vacuum, 

A(k 2 ) - 



-1 



k 2 — m 2 + ie 



(A.« 



(A.9) 



(A.10) 



The thermal propagator may be diagonalised in the form 

D ab (k ,k) = C/ ac (/%o)[diag{A(/% ,fc),-A*(A : o,fc)}] c ' i C/ db (A : o) 
with the elements of the diagonalising matrix as 

U n = U 22 = VTT^, u 12 = u 21 = V^ 

Writing spectral representations, one can show that U diagonalises not only the free propagator, but also the 
complete one [12j [18] . 

Free vector propagator 
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The thermal propagator for a massive, spin-one particle may be derived in a similar way. Denoting its field 
by Pfj,(x), it has the propagator 

D^(t,t';x -x' ) = i(T c p^(T,x)p u (T',x')) . (A. 11) 

satisfying the differential equation 

[-^(□ c + m 2 ) + d^d v ]D v P{ry ] x-x') = <5 c (r - r')5\x - x') (A.12) 

The thermal boundary condition is again given by an equation similar to that of Eq.(A.4). As before we take 
the spatial Fourier transform, when the the tensor components satisfy 

u 2 D 00 (t - t', k) + iki-^-D l0 (T - t', k) = -5{t - t') 
(it 

uj 2 D°Ht - t', k) + iki—D l Hr - t', k) = 
ar 



9,1 ( ^2 +uj2 ) + k ' k > 



D l Hr - t', k) - iki^D°i = g>5{ T - r') (A.13) 
ar 



It is easy to check that these equations as well as the boundary condition is solved by 



Dij(T 


-r',k) = 


f kikj ^ 
m 2 J 


D(t- 


D 0i (T 


-r',k) = 


.k d 

l—^-D{T-T 

m z ar 


,k) 


D 00 (t 


-r',k) = 


(l+ 1 

\ m 2 dr 2 


)d(t 



\k) (A.14) 

where D satisfies the same Eq. ()A,6P as for the scalar field. We now take the temporal Fourier transform as 
before and assemble the resulting four dimensional Fourier Transform of the four components of the thermal 
vector propagator in terms of that for the scalar propagator 

Df v (k) = f-g^ + D ab (k) (A.15) 

the Lorentz structure remaining the same as for the vacuum propagator. 

Diagonalisation of complete vector propagator 

We now turn to Dyson equation (3.3) in the text. As already stated above, both the free and the complete 
propagator can be diagonalised by the same thermal matrix U ; thus 

Gf u (k , k) = U ac (k )[di a g{G^(k , k), -%,{ko, k)}] cd U db (k ) (A.16) 

As a consequence, the matrix II?* is also diagonalisable by (U~ 1 ) ab , 

(9) = [U-Hqo)] ac [^g{Ti, u (q), -%M}] cd [U-\ qo )] db (A.17) 

We thus get the Dyson equation (3.4) in the text for the barred quantities, that are free from the thermal 
indices. The diagonalised tensor can be related to any of the components of the corresponding matrix. For 
example, it is related to the 11 component as 

ReV = ReU jl 

ImTT^ = tanh(/3w/2)Imn^ (A.18) 
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